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We investigate the low transverse momentum heavy quarkonium (J/ip, T) productions in the 
forward rapidity region of pp and pA collisions at the LHC as an important probe to the transverse 
momentum tomography of the gluons in hadrons in the Color Glass Condensate (CGC) framework. 
By implementing the Sudakov resummation consistently in the CGC formalism, we achieve an 
excellent agreement between the improved CGC calculations and the LHC data. We show that both 
the small-x and the Sudakov effects are essential for a complete description of heavy quarkonium 
productions in the low transverse momentum region. This provides a solid foundation to study the 
small-a; gluon saturation in a big nucleus with the future pA programs at the LHC. 
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Introduction- Due to the rich and complex dynamics of strong interaction and its close connection to various 
interesting phenomena, heavy quarkonium production in high energy collisions has been one of the most intriguing 
subjects in the study of strong interaction physics. In particular, as the objective of this paper, it has been shown 
that the low-pj_ spectrum of forward rapidity heavy quarkoniums produced in high energy pp (pA) collisions carries 
important information of the gluon saturation at small-a;. The concept of gluon saturation |l|, 2] comes from the 
anticipation of gluon recombination when its density becomes extremely high in high energy collisions. The associated 
dynamics is one of the challenging physics topics in current and future nuclear science program around the world Q. 

There have been intensive theoretical studies on heavy quarkonium productions in the small-a; saturation framework, 
i.e., the Color-Glass-Condensate (CGC) formalism. They have shown that the typical transverse momentum of 
produced quarkonium is around the saturation momentum Q s {x g ) EHH . As a function of the longitudinal momentum 
fraction of target hadron x g , the saturation momentum increases when x g decreases, and it characterizes the typical 
transverse momentum that small -x gluons carry in a hadron. It can also be derived from the boundary which separates 
the dilute regime from the dense saturated regime where gluons evolve non-linearly. The non-linear evolution equation 
at small-x (or high energy) is known as the Balitsky-Kovchegov (BK) equation 13, 14], which has been implemented 
numerically in CGC formalism. 

The physical picture of quarkonium productions in the CGC framework is as follows. Either before or after 
interacting with the dense gluons in the target hadron, a large x gluon from projectile proton fluctuates into a pair 
of heavy quarks, which eventually forms a bound state of quarkonium with transverse momentum p±. In the low p± 
region, the main contribution of the transverse momentum broadening is due to gluon saturation effects. At forward 
rapidity y , the kinematics indicates that x g from the target side is much smaller than the one in the projectile. This 
is the ideal region for the study of small-cc physics. Therefore, the heavy quarkonium productions in the forward pp 
and pA collisions can provide an important probe for the gluon saturation at small-a:, whose tomography imaging has 
attracted strong interests in hadron physics community [3]. 

Thanks to the abundant data ]15l - |23j and wide kinematical coverage at the LHC, we are able to test the 
CGC/saturation physics calculations for pp and pA collisions by comparing them to experimental data. Early 
studies of J/A productions [§1 ITol - [T2| have shown that the CGC can provide decent description of the data, especially 
the transverse momentum distribution of the produced J/A. Due to similar kinematics and physical parameters 
except for differences in mass, one expect that the p± broadening of T should be similar or slightly smaller than 
the one observed for J/A due to saturation effects [l2|. However, the LHCb experiments |la Oil |2Q[ have found 
consistently larger mean p± for T for a wide range of rapidity windows, as compared to that of J/tp, as well as the 
theoretical calculations of [l2|. In this paper, we will solve this puzzle. 

We will show that one needs to take into account the Sudakov resummation together with the small -x effect, 
in order to consistently describe the J/tf and T data. In a recent paper [24j, it has been demonstrated that the 
Sudakov resummation can be performed simultaneously with the small-x resummation. Furthermore, Refs. f25l |26| 
computed the corresponding Sudakov double logarithms (in terms of ol s Ca In ' 2 (M 2 /p 2 j_) with M the quarkonium mass) 
for heavy quarkonium productions in hadron-hadron collisions in the collincar/CGC hybrid formalism. The results 
are consistent with those obtained in the Collins-Soper-Sterman (CSS) formalism j27j. Based on the commonly used 
CGC calculations for heavy quarkoniums, we implement the Sudakov effects and compute the p± spectrum for both 
J/tf and T. Since the typical transverse momentum of the produced quarkonium, namely the saturation momentum 
Q s ~ 1-2 GeV, is not much smaller than the J/tf mass, the Sudakov effect is found to be negligible for J/if> productions 
at the LHC. On the other hand, the Sudakov logarithms are important for T which consists of bottom quark with 
the mass around 4 GeV. By performing the numerical calculation with the Sudakov factors, which provide extra p± 
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broadening for T spectra, we find a complete agreement with the LHC data for both J/ip and T in the low-px region. 
This shall provide an important step forward to study the small-a: gluon saturation in the future pA program at the 
LHC. 

The implementation of the Sudakov resummation in the CGC framework- In the forward rapidity region, transverse 
momentum of gluon coming from incident proton should be of order less than O(Aqcd), which allow us to apply the 
usual collinear gluon distribution for the dilute incident proton side. Therefore we can start with the heavy quark 
pair production cross section formula in pp (pA) collisions for minimum bias events at leading order [§j as follows 
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where q q ± (q q ±) and y q (y q ) are transverse momentum and rapidity for produced quark (antiquark) respectively. 
x\ (x 2 ) is the longitudinal momentum fraction of the dilute proton (the dense nucleus) carried by the incoming 
gluon. For quarkonium production, we find x\, 2 = \/M 2 + p\/yfse within 2 — > 1 kinematics where M is invariant 
mass of a quark and antiquark pair and p± = q q _ j_ + q q ± is transverse momentum of the pair. p± is nothing but 
quarkonium transverse momentum. Y = l/21n((g+ + q^)/(qf + q/f)) is the pair rapidity. X\G is the usual collinear 
gluon distribution function with p being the factorization scale. In the above formula, we have k 2 ± = p± due to 
momentum conservation. The hard matrix element is given by 3 co u = S®®(®® + S®®j® + 3®’^ with 
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where a± = q q j_ — k±. S®®(®® represents the scattering of heavy quark pair produced from gluon splitting off the dense 
nucleus and S®^ corresponds to the scattering of gluon before heavy quark pair splitting off the dense nucleus, whereas 
5®®jf is derived from the interference between the above two cases. In the above expressions, we have introduced 
the momentum fraction z for the quark z = qff / (g+ + q/f ) and z q = 1 — z for the antiquark. The multipoint point 
function <p A y constructed with the color singlet dipole amplitude as follows; 


^ g (k 2± ,k ± ) 


= S i 
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: S± Fy g (k±) F Yg (k 2 ± — k_ l), 


( 3 ) 


where we denote S± = S± Nc ^ ± with S± = ttR\ as the transverse area of target nucleus and Fy (k±) = 

f ( 2 ttF e ~ lk± x± Sy 3 (x± ) . Here Y g = ln|^ is the evolution rapidity of the gluon in the target with x 0 = 0.01. 
The rapidity evolution of Sy is described by the nonlinear BK equation. For numerical calculations, we use the 
BK equation with running coupling correction (rcBK) [2|j. The initial condition of the rcBK is the so-called 
McLerran-Venugopalan model [29j with modified anomalous dimension at xq = 0.01: 


SV^oO-l) = exp 




( 4 ) 


It is well known that HERA e+p global data fitting foot - fri provides a nice constrained initial condition by using 
Eq. ([4]) . In this paper, since the results are not sensitive to the initial condition, we adopt the same parameters set 
for Eq. (@| as found in Ref. Q. In pA collisions, we change the initial saturation scale as Q 2 A 0 = cA 1 / z Q 2 0 with 
c = O.5 0. 

The differential cross section of the pair production with the Sudakov factor in momentum space is given by 
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where we denote the Fourier transform of the Sudakov factor with the gluon pdf as follows 

F1 


Fsud(M,l ± ) = j 


( 7 ) 


where eg = 2e 7E with 7 e the Euler-constant and the factorization scale is chosen as p = 
In Eq. dTJ) , the exponent of the Sudakov factor is given by 


<Ssud(Af, b) = S , perp(M, 6*) + Snp(M, b ) (8) 

where we have introduced b+ = b/\J 1 + (6/6 max ) 2 by following the CSS formalism [27j to separate the perturbative 
part (b < femax) from the nonpertubative part (b > 6 max ). 6 max is a cutoff scale to separate these two regions. 
SperpiQ, b) is calculable perturbatively at small b which can be cast into 


S peip (M,b) 




(9) 
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respectively. For the one loop correction, they are given by A W = Ca and B W = — {bg + \Ssc)N c where bg = 
N c — -A-. As for B^\ the factor <5 8c is significant only in the production of a color octet quark-antiquark pair. 

In our calculation, because the quarkonium is produced by the color octet channel in the color evaporation model 
(CEM) with large-A^ c , we will take into account the contribution from this term. 

For the non-perturbative part Snp(/i, 5), we adopt the one from Ref. fildl derived for the color octet quarkonium 
production as follows 


Snp ( M , b) = exp 


Y (-Si -S2‘n (^) - 


gig 3 ln( 10 Ctei:E 2 ) 


( 10 ) 


where g\ = 0.03, g 3 = 0.87, and <?i x g 3 = —0.17 are obtained by the data fitting within the NRQCD factorization 
with Q o = 1.6 GeV and 6 max = 0.5 GeV is originally fixed in Ref. [34j. Since we only employ the collinear gluon 


distribution for the large x gluon from the incident proton, Eq. (flOl) differs by a half as compared to Ref. 34]. As found 
in Ref. [HJ , all the parameters are determined by assuming that the color octet heavy quark pair is dominant channel 
of the quarkonium production at low transverse momentum. In the same spirit, we can use the same parameters for 
the CEM. 


Numerical results- Now we turn to the numerical calculations for quarkonium production in the forward pp and pA 
collisions by adopting the CEM. We expect that the Sudakov effect should also apply in other models. In the CEM, 
the invariant mass of heavy quark pair is integrated from the bare quark pair mass to the threshold of decay into 
an open heavy flavor meson pair, namely 2m < M < 2 Mh- We fix m = 1.2 GeV and Mh=D = 1.864 GeV for J/^> 
production, and m. = 4.5 GeV and Mh=B = 5.280 GeV for T respectively. The produced quark pair is going to be 
bound into a quarkonium with the probability F q q^Q. This empirical factor is interpreted as a normalization factor 
because it depends on quark mass and factorization scale and etc. We have also included any K factor in association 
with higher order correction in F qq ^Q as well. 

Some setup parameters in numerical calculations are listed here. We fix R p = 0.9 fm and Ra = 8.5 fm, a s = 0.3 
in the hard part. As for the running coupling a s (Q) in S peTp (M, 6*), we use Aqq D = 0.173 GeV for 4-flavor in T 

/ o\ 

productions and Aq ^ = 0.214 GeV for 3-flavor in J/tf) productions. These values are extracted to reproduce the 
low-Q data of a s (Q) (3]| up to Q = 10 GeV by using 1-loop running coupling. In addition, we adopt CTEQ6 pdf Q 
in numerical calculations: CTEQ6L for Eq. JT|) and CTEQ6M for Eq. d^]). 

The main focus of this paper is on the transverse momentum broadening of T production, while the parallel 
calculation for J/ij) production is mainly treated as a comparison to illustrate the Sudakov effect. We show in FIG.|T| 
the differential cross section of prompt J/^/i production as a function of p± in pp collisions at y/s = 7 TeV and at 
forward rapidity. The results of Eq. m, which is without the Sudakov effect, are computed with only small -x evolution 
by using rcBK evolution. These results are roughly in agreement with the data, especially for Y = 4.25, where the 
hybrid formalism is supposed to work the best. With the Sudakov factor, we find that Eq. ([5]) provides slightly better 
agreement as compared with the result of Eq. ©. In general, the additional p± broadening provided by the Sudakov 
factor is negligible especially at Y = 4.25 where the saturation effect is the strongest. We would like to emphasize 
that the Sudakov resummation is not applicable when p± is larger than the hard scale of the system. When p± ~ AI, 
we should switch to the fixed order CGC calculation, which is responsible for the large p± region of the spectrum. 
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FIG. 1. Double differential cross section of J/'ij) as a function of p± for Y = 2.25, 3.25, and 4.25 in pp collisions at fs = 7 
TeV. Blue solid line is obtained by using Eq. m and the uncertainty band is coming from a change of factorization scales 
(2 < p < 30 GeV) for the collinear gluon distribution function. Red solid line denotes the result of Eq. [5] at 6 max = 0.5. We 
choose Fj/^ = 0.0975 for Eq. JlJ and 0.1495 for Eq. ([5]). The LHCb data for prompt production is taken from Ref. ficj] . 
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FIG. 2. Double differential cross section of T(1S) multiplied by a branching ratio of Y(1S) decay into a lepton pair as a 
function of p± for Y = 2.25, 3.25, and 4.25 in pp collisions at fs = 7 TeV (solid lines) and 14 TeV (dotted lines). We choose 
Fr(is) — 0.488 for Eq. JIJ and 0.390 for Eq. (0. The LHCb data is taken from Ref. fill . 


FIG. [2] shows the results for T production at fs = 7 TeV and 14 TeV. As shown in Refs. @, HH, the T production 
with only the small-a; resummation does not yield enough transverse momentum broadening as compared to data. 
We can find that the results from Eq. ([5]) reproduce the data points very well due to the gluon cascade characterized 
by the Sudakov factor. The peak of the spectrum from Eq. m is located around p± = 1 GeV while the Sudakov 
factor shifts the peak position to around 3 GeV. Of particular importance is that the Sudakov factor in association 
with large M gives additional strong broadening of the p± distributions for T production. The Sudakov resumma¬ 
tion is indispensable in order to consistently describe all heavy quarkonium productions within the CGC/saturation 
formalism. 

In FIG. El we finally show the differential cross section of T production as a function of p± in pPb collisions at 
yfs = 5.02 TeV and fs = 8.8 TeV by setting Q 2 A 0 = 3 Q 2 s0 in the initial condition. We find that the Sudakov effect 





p ± (GeV/c) p ± (GeV/c) P ± (GeV/c) 

FIG. 3. Double differential cross section of T(IS') in pPb collisions as a function of p± for Y = 2.75 (solid lines) and 3.75 
(dotted lines) at fs = 5.02 TeV (Left) and fs = 8.8 TeV (Middle). Right: Nuclear modification factor of Y(IS') as a function 
of p± at 5.02 TeV and 8.8 TeV for Y = 2.75 (solid lines) and 3.75 (dotted lines). 
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in pA collisions is less pronounced as compared to pp collisions, since the saturation effects becomes much stronger 
in pA collisions due to the nuclear enhancement and it starts to dominate over the Sudakov effect. We also show 
in FIG. [3] the T nuclear modification factor R p pb = ^ pPh as a function of p±. The deviation of R p pb from unity 
can reveal the presence and strength of the gluon saturation effect. As shown in the R p pb plot, the suppression gets 
stronger in higher energy and more forward rapidity due to increased saturation effect. 

Conclusion- In this paper, by computing the p± spectra for forward J/'F and T productions at the LHC, we 
demonstrate that the Sudakov resummation on top of the small-x evolution in the saturation formalism is essential to 
interpret the LHC data. Further theoretical and experimental efforts along this line can help us quantitatively study 
the k± dependent gluon distributions inside high energy protons and nuclei, and provide us insightful information 
about the hadron tomography especially in the low-a; limit, as well as important evidences for the onset of gluon 
saturation |3j. 
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